† Corresponding author. E-mail:
Dynamics of the Au + H2 reaction are studied using time-dependent wave packet (TDWP) and quasi-classical trajectory (QCT) methods based on a new potential energy surface [Int. J. Quantum Chem.
Since the excellent catalytic effect of gold nano-cluster on low temperature oxidation reaction of CO was found by Haruta and co-workers[1] in 1989, a number of experiments have been performed on nanoscale gold clusters. In addition, properties of nanoscale gold clusters such as electronic, optical, chemical, and catalytic were reported in previous theoretical studies.[2–9] To further explore the properties of gold and its clusters, a lot of theoretical studies have been carried out for the reaction of transition metallic element Au and its ions with H2 molecules.[10–20] Using the complete active space self-consistent field (CASSCF) and configuration interaction method, two low-lying electronic states of AuH2 systems were studied by Balasubramania and Liao[10] in 1988. The group of Balasubramanian[11] restudied the AuH2 system in 2004. In their work, the bending potential energy surface (PES) was studied using the CASSCF/MRSDCI methods on three low-lying electronic states of AuH2 systems and they found that there is a barrier on the path when the ground state AuH2 system decompose into Au(2S) + H2 channel. Using highly correlated wave functions with atomic pseudopotentials and spin–orbit interactions, the electronic structure, reactivity, and spectroscopy of the AuH2 system were investigated by Chambaud et al.[12] in 2005.
As far as we know, there are few studies on dynamics of the Au + H2 reaction due to no global PES of AuH2 systems until 2010. The first global PES of the ground state of AuH2 system was constructed by Zanchet et al., called the Zanchet PES.[13] In the ab initio calculation, multi-reference configuration interaction (MRCI) method with one-electron Gaussian-type basis set ECP60MDF[21] and augmented correlation-consistent polarization valence triple zeta (AVTZ)[22] basis set were employed. Then, the PES was constructed through GFIT3 C procedure[23–25] by fitting 2376 ab initio points. Based on the Zanchet PES, some dynamics studies have been performed. For example, the reactions of Au with H2 and its isotopic variant HD, D2 were studied by Yuan and co-workers[14,15] using the time-dependent wave packet (TDWP) and quasi-classical trajectory (QCT) methods. The dynamic properties such as reaction probability, integral cross section (ICS), and differential cross section (DCS) were studied at state-to-state level of theory. The Zanchet PES with relative little ab initio points and the basis sets used in the ab initio calculation also have a space to improve. Thus, Yuan and co-workers rebuilt the PES of the ground state of the AuH2 system by fitting 22853 ab initio points with neural network method (YLYC PES).[16] In their work, the dynamics calculation was also performed at the initial state (v0 = 0, j0 = 0) and the reaction probability, ICSs and DCSs were calculated in the collision energy range from 0.01 eV to 1.0 eV. In addition, there are also some studies on the reactions of gold ions with H2 molecules.[17–20]
As discussed above, the reaction of Au with H2 and its isotopic variant HD, D2 were well studied based on the Zanchet PES. However, there are few theoretical studies of the Au + H2/HD/D2 reaction based on the new YLYC PES which is more accurate than the Zanchet PES at state-to-state level of theory. Thus, the aim of the present work is to perform the dynamics calculation at state-to-state level of theory and to further understand the reaction mechanism. This paper is organized as follows: the theory of TDWP and QCT methods are presented in Section
As one of powerful computational tool of chemical reaction dynamics, the TDWP method is widely used in the atom-diatom reaction[26–33] and the reaction involving polyatomic molecules.[34,35] For the Au + H2 reaction in the body-fixed reaction Jacobi coordinate, the Hamiltonian for a given total angular momentum J can be written as
Before the dynamics calculation, an initial wave packet with specific initial state should be set up. The initial wave packet consists of three parts: a Gaussian wave packet in the space-fixed frame, a ro-vibrational eigenfunction of the H2 molecule, and a total angular momentum eigenfunction with parity of system ε = (−1)j0 + l0. The initial wave packet is presented as
The reactant coordinates based method[37,38] is used to extract the state-to-state S-matrix information. The body-fixed time-independent scattering wave function can be obtained through an orthogonal transformation matrix
Then, through a standard transformation, the state-to-state scattering matrix is transformed into the helicity representation
Based on the YLYC PES, standard QCT[39] calculations for the Au + H2 (v0 = 0, j0 = 0) reaction were carried out in the collision energy range from 1.5 eV to 3.0 eV, with steps of 0.1 eV. The trajectories started at a distance between the incoming atom and the center of mass of the H2 molecule of 8.0 Å, and 1000000 trajectories are sampled for each collision energy. In the calculation of trajectories, the integration step is set up as 0.02 fs to ensure the numerical stability. For the total reaction probability of total angular momentum J = 0, the impact parameter (b) is set as 0 and the reaction probability can be obtained by the following formula:
The YLYC PES is used in the present TDWP and QCT calculations. The properties of the YLYC PES are listed in Fig.
The parameters used in the TDWP calculation have large effect on the final results, thus the numerical parameters should be tested before calculation. A number of tests were performed on the total reaction probabilities of total angular momentum J = 0. Finally, the converged parameters were obtained and also used in the J > 0 calculation. The parameters are presented in Table
The total reaction probabilities of the title reaction at the total angular momentum J = 0 are listed in Fig.
The total and vibrational state-resolved ICSs are presented in the Fig.
The ro-vibrational state-resolved ICSs of the title reaction calculated by the TDWP method are presented in Fig.
The total DCSs of the Au + H2 reaction at four selected collision energies 1.8 eV, 2.2 eV, 2.6 eV, and 3.0 eV are plotted in Fig.
In order to further understand the reaction mechanism, the vibrational state-resolved DCSs of the title reaction at the collision energies of 1.8 eV, 2.2 eV, 2.6 eV, and 3.0 eV are displayed in Fig.
In order to further understand the reaction mechanism, the rotational state-resolved DCSs of the vibrational state v′ = 0 at the four selected collision energies. As seen in the figure, at low collision energy of 1.8 eV, it is almost forward-backward symmetry scattering. However, the side-way signals become more and more apparent when the collision energy increases and the peaks are around 90° at a relative high rotational quantum number.
The initial state v0 = 0, j0 = 0 dynamics calculation of the Au + H2 reaction was carried out in the collision energy range from 0.01 eV to 1.0 eV. The TDWP method with the second order split operator and the QCT method were used in the dynamics calculations. The reaction probability, ICS, DCS, the ro-vibrational state-resolved DCS and ICS were reported and compared with previous theoretical studies. Some discrepancies can be found between previous values and the present results due to the different PESs used in the calculations. Furthermore, three reaction mechanisms were found in the dynamics calculations, the complex-forming reaction mechanism is dominated in the low collision energy region, and the abstract reaction mechanism plays a dominant at the high collision energies. In addition, the side-way scattering signals were also found when the collision energy further increases.
[1] | |
[2] | |
[3] | |
[4] | |
[5] | |
[6] | |
[7] | |
[8] | |
[9] | |
[10] | |
[11] | |
[12] | |
[13] | |
[14] | |
[15] | |
[16] | |
[17] | |
[18] | |
[19] | |
[20] | |
[21] | |
[22] | |
[23] | |
[24] | |
[25] | |
[26] | |
[27] | |
[28] | |
[29] | |
[30] | |
[31] | |
[32] | |
[33] | |
[34] | |
[35] | |
[36] | |
[37] | |
[38] | |
[39] |